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Collisions between cosmic bubbles of different vacua are a generic feature of false vacuum eternal 
inflation scenarios. While previous studies have focused on the consequences of a single collision 
event in an observer's past, we begin here an investigation of the more general scenario allowing for 
many "mild" collisions intersecting our past light cone (and one another). We discuss the general 
features of multiple collision scenarios and consider their impact on the cosmic microwave back- 
ground (CMB) temperature power spectrum, treating the collisions perturbatively. In a large class 
of models, one can approximate a multiple collision scenario as a superposition of individual colli- 
sion events governed by nearly isotropic and scale-invariant distributions, most appearing to take 
up less than half of the sky. In this case, the shape of the expected CMB temperature spectrum 
maintains statistical isotropy and typically features a dramatic increase in power in the low multi- 
poles relative to that of the best-fit ACDM model, which is in tension particularly with the observed 
quadrupole. We argue that this predicted spectrum is largely model-independent and can be used to 
outline features of the underlying statistical distributions of colliding bubbles consistent with CMB 
temperature measurements. 



I. INTRODUCTION 

The theory of cosmological inflation has been very suc- 
cessful in proving potential answers to several deep ques- 
tions in early universe cosmology, as well as testable (and 
tested) cosmological predictions. However, a crucial "side 
effect" of many versions of inflation is that the exponen- 
tial expansion tends to continue on eternally into the fu- 
ture, with only pockets of the universe ceasing to inflate 
locally. This scenario is known as eternal inflation (see 
e.g. Refs PU3] for reviews). 

Eternal inflation can be driven by several mechanisms. 
In the present study, we concern ourselves with the "false 
vacuum" variety, in which inflation occurs as a result of 
an inflaton (usually taken to be an effective scalar field, 
f) being trapped in a metastable false vacuum. Infla- 
tion ends locally when the field transitions to a nearby 
"truer" vacuum on the potential. If the probability per 
unit four- volume of such a transition, A, is small com- 
pared to the Hubble scale of the inflating false vacuum 
Hp (as typically assumed), then false vacuum inflation 
is eternal. For many potentials, the dominant transition 
mechanism is bubble nucleation [3] , which can quite nat- 
urally produce open FRW universes much like our own. 
In this study, we will work under the assumption that 
we live in one such universe, dubbed "the observation 
bubble". 

A fascinating consequence of false vacuum eternal in- 
flation is that other bubbles inevitably form in the four- 
volume to the past of the observation bubble wall. These 
bubbles collide with the observation bubble wall, poten- 
tially leading to detectable signatures of "cosmic bub- 
bles" outside of our own, provided that an observer can 
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exist to the future of the collision. Much recent work 
has been dedicated to the study of such relics [5H22] (for 
a general review, see e.g. Ref. [23] )■ In particular, 
it has been shown that if a collision is "mild" enough, 
it can potentially be compatible with the formation and 
survival of observers in the post-collision region [7J [8] 
while producing testable predictions for observations of 
the Cosmic Microwave Background (CMB)[Tfjl IT71 [TTj]. 
bulk galactic flow [12], and possibly other cosmological 
observables (see e.g. Ref. [23] for a discussion). 

While the observation bubble wall undergoes a diver- 
gent number of collisions, whether or not these collisions 
lie within the observer's past lightcone at decoupling de- 
pends on several factors, such as the position of the ob- 
server and the cosmology inside the observation bubble, 
which we review in more detail in Sec. [Hi Previous stud- 
ies have focused on the outcome of a single observable 
collision, or a few collisions whose affected post-collision 
regions do not intersect. In these cases, the azimuthal 
symmetry of the collisions' perturbations to the obser- 
vation bubble spacetime is maintained and one expects 
to see an azimuthally symmetric disk on the CMB sky 
corresponding to each collision [10l [17] . Recent studies 
[TO1 [50] have searched for these disk-like patterns in the 
WMAP7 CMB temperature data and have found a few 
candidates for collision events. 

As we will see in Sec. [IT] however, for arbitrary choices 
of the various parameters, the number of observable bub- 
bles N can potentially range from zero to very large num- 
bers, and there is no particular reason to expect ~ one 
observable collision rather than none or many. When 
N ^> 1, the picture differs from that of previous studies 
in several respects: 

1 . The future light cones of the different collisions will 
tend to intersect one another. Thus, the effects aris- 
ing from various collisions can potentially interact. 
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2. One must address the possibility that some bubbles 
may in fact collide with each other before impacting 
the observation bubble wall, leading to a qualita- 
tively different collision scenario than before. 

3. Since any given part of an observer's LSS will have 
been affected by several different collisions, the az- 
imuthal symmetry of an isolated collision is no 
longer evident. 

4. With a large number of bubble collisions, one may 
be able to make meaningful statistical predictions 
about what a given observer can expect to see. 

The purpose of this study is to explore these issues and 
take the first step towards assessing how many-bubble 
scenarios can be constrained by current and future ob- 
servations. As a first look, we will focus primarily on 
the effect of the collisions on the CMB temperature, de- 
ferring an analysis of the other relevant observables to 
future study. 

Clearly if we are to be living in such a large- N sce- 
nario, the effect of each collision must in some sense be 
"small" so as not to disrupt the subsequent cosmological 
evolution in the observation bubble. As we will discuss 
in Sec. jlVj the strength of a collision is determined by 
the position of the nucleation center of the colliding bub- 
ble in the exterior de Sitter (dS) space, as well as the 
details of the underlying potential, which are unknown. 
Thus, in analyzing scenarios with N > 1, we will simply 
start with the assumption that an observer exists with 
a large number of collisions intersecting the visible por- 
tion of the last scattering surface and study what might 
be observed. This is in keeping with previous work on 
the CMB temperature profile for single collisions. We 
do not consider scenarios in which some collisions might 
prevent inflation to the future of the collision, such as 
those in small-field type potentials (see the discussion in 
Ref. [23]). When discussing an observer's expectations, 
we therefore mean the expected value considering only 
mild collision scenarios. 

The current study is organized as follows: in Sees. \H[ 
|IV| we set up the various parts of the problem and argue 
that many scenarios with multiple bubble collisions can 
be considered as a superposition of individual impacts, 
each appearing to take up less than half of the observer's 
sky. In Sec. [V] we compute the expected signal in the 
CMB temperature power spectrum. Sec. |VI| contains our 
conclusions. 



II. BUBBLES, COLLISIONS, AND OBSERVERS 
A. The Setup 




FIG. 1. Conformal slice of the basic bubble collision setup 
in the so-called "steady state frame" . The observation bub- 
ble shown is pure dS with Hf/Hi ^> 1. The region out of 
which colliding bubbles is labelled as such, and the red line 
demarcates the initial value surface. 



coordinates X^ and metric 77^, where /j, v = 0, ...,4. 
The most useful coordinatization of dS for our purposes 
will be the so-called "global slicing", with coordinates 
(T, 77, 9, <f>) defined in the embedding space by 



X Q = HpHanT 

F cosT 
Xa = Hp 1 



(1) 



cos 77 



F cosT' 

where — 7r/2 < T < ir/2, < 77 < tt, and u)i 
(cos#,sin#cos</>,sin$sin0). This induces the metric 



ds 2 = 



1 



H 2 cos 2 T 



(-dT 2 + drj 2 + sin 2 77 dfl 2 ,) 



(2) 



which covers the entire dS manifold. The 4-volumc cle- 
ment in these coordinates is 



dV A = H 



_ 4 smi| 



cos 4 T 



dTdr) dQ,2- 



(3) 



The basic setup we consider is the same as that of 
Refs. 6-8, 23 , to which we refer the reader for further 
details. A small, thin-walled observation bubble is nucle- 
ated out of the false vacuum via the Coleman-deLuccia 
(CDL) process |lj and subsequently expands with a tra- 
jectory which we approximate as null{^| The observation 



We begin with a false vacuum, which we take to be dS 
space with associated Hubble constant Hp. We can de- 
scribe the false vacuum dS as a hyperboloid r\ liV X li X 1 ' = 
Hp 2 embedded in a 5-dimensional Minkowski space with 



We note that, while this small-bubble, thin-wall limit S\ is conve- 
nient and used here for both collision and observation bubbes, it 
is not necessarily natural and future work is needed to go beyond 
this approximation. 
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bubble starts off curvature dominated and undergoes a 
period of subsequent slow-roll inflation, which we model 
as dS space with Hubble constant Hj. Throughout this 
study, we assume Hp /Hi 3> 1 except where noted; this 
seems natural, as Hi < Hp for tunneling to occur, and 
there is no obvious reason why the scale of inflation in the 
observation bubble should be similar to the scale associ- 
ated with the false vacuum. After inflation, we assume a 
standard ACDM cosmology, although we typically ignore 
matter domination and the possibility of late-time accel- 
erated expansion in our calculations, which should induce 
small corrections. This scenario is sketched in Fig.[T]for a 
conformal slice of 8, <fi =const. (Note that as drawn, the 
observation bubble is pure dS with Hp /Hi 3> 1). In this 
diagram the bottom (top) boundary of the false vacuum 
corresponds to T = —ir/2 (tt/2), while ij ranges from 
to 7r from left to right. 

The observation bubble interior spacetime is given 
by the analytic continuation of the CDL instanton [3], 
which, in our approximation, yields an open FRW cos- 
mology inside the null cone emanating from the nucle- 
ation center, with metric 



ds 2 



-dr 2 + a 2 (r) (d£ 2 + sinh 2 



(4) 



Spacelike surfaces of constant r are 3-hyperboloids, Ha, 
reflecting the SO(3, l)-invariance of the CDL instanton. 
The observation bubble and false vacuum dS spacetimes 
can be joined together by gluing across the bubble wall, 
as described in Ref. [Hj. 

Along with the observation bubble, other bubbles can 
also nucleate out of the false vacuum with associated 
probability per unit 4- volume Aj, where i represents the 
bubble's corresponding vacuum. The region of dS that 
can nucleate colliding bubbles is labelled in Fig. [T] in the 
"steady state" frame in which the distribution of bub- 
bles is statistically independent of position and time, and 
the conformal diagram is independent of angle. Bubbles 
nucleated above the top boundary of this region never 
intersect the observation bubble wall, while the lower 
boundary ensures that the observation bubble was in- 
deed nucleated out of the false vacuum and not some 
other bubbleFJ Bubbles formed to the right of the red 
boundary will have future light cones that encompass fu- 
ture infinity of the false vacuum dS space, and so as in 
previous studies we postulate a "no bubble" initial con- 
dition surface as shown. 

While the diagram in Fig. [T] depicts the situation in 
the steady-state frame, there are other reference frames 
that will be useful to consider. One can define a frame 
in which the observer is at the center of the observation 
bubble £ = (the "observation frame" ) , as well different 
"collision frames" in which the observation bubble and a 



given colliding bubble nucleate at the same global slic- 
ing time T = 0. To boost between these different frames 
in both the observation bubble and exterior dS, one can 
apply the appropriate Lorentz transformation in the em- 
bedding space PJ As discussed in detail in Ref. [H] , in the 
observation frame the distance from the observation bub- 
ble wall to the initial value surface is angle-dependent. 
This fact plays an important role when computing the 
distribution of bubbles an observer expects to see, as we 
will discuss in the following subsection. 

Since a collision preserves an SO (2, 1) symmetry, it can 
be fully described in 1+1 dimensions, with each point 
representing a 2-hyperboloid with line element dHi . The 
dS false vacuum can be conveniently foliated according 
to this symmetry, yielding the metric 



ds 2 



(1 + H" l z 2 Y x dz 2 + (1 + H 2 z 2 )dx 2 + z 2 dH\. (5) 



We can then place the observation and colliding bub- 
ble nucleation centers along the rr-axis, and describe the 
kinematics of the collision by (x c , z c ), the nucleation cen- 
ter of the colliding bubble (the observation bubble is at 
the origin). In the collision frame, this simplifies fur- 
ther: since both bubbles, by construction, nucleate at 
z = 0, all of the kinematics are described by z c , which 
corresponds uniquely to the de Sitter-invariant distance 
(rjc in the collision frame) separating the bubbles. This 
quantity is responsible for determining the strength of a 



given collision, as we will discuss in Sec. IV 



B. Probabilities And Observer Expectations 

Ascertaining the general features of a multiple colli- 
sion scenario requires information about the various pa- 
rameters that enter into the theory - in particular, the 
expected number of observable bubbles intersecting the 
observer's last scattering surface, as well as the distribu- 
tions of the angular sizes and strengths of the collisions 
that would enter into e.g. the CMB temperature profile. 
This section reviews the extensive recent progress made 
in understanding these distributions, focusing on what 
is relevant to the task at hand and offering some novel 
insights. We refer the reader to Ref. [53] for a good sum- 
mary and for details of many of the calculations discussed 
here. 

We begin in the observation frame, in which £o = 0. 
Under our assumptions, the number of bubbles intersect- 
ing part or all of the observer's last scattering surface is 



N = \V, 



obs 



(6) 



where V£ is the 4- volume in the false vacuum dS space- 
time out of which such a bubble can nucleate, computed 



2 Throughout this study we work in the approximation in which no 
bubbles can nucleate within another bubble. This approximation 
is discussed in more detail in Sec. IIVI 



3 One should note that in general once there are collisions it will 
take more dimensions to embed the spacetime and boost. 
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by integrating the 4- volume element Eq. [3] over the rele- 
vant region. There are two regions in the exterior dS in 
which (IV4 can potentially greatly exceed its natural scale 
of Hp A , and hence contribute the most to V^ obs and the 
total count of bubbles. 

First, for observers with £0 — ► 00 j there is a large con- 
tribution to N from portions of V^ 01 " 3 near past infinity 
of the false vacuum. This results in an anisotropic peak 
in the distribution of bubbles in the direction towards 
the initial value surface (defined as 9 = 0), which can 
be thought of as the time dilation of the observer's time 
relative to the "steady state" cosmological time. This 
anisotropy was first recognized in Ref. [5] and dubbed 
the "persistence of memory" . Bubbles nucleated in the 
region near past infinity will cover essentially the full sky 
and enter the observer's past light cone at small r and 
were thus dubbed "early time" collisions in Ref. [6]. As 
discussed in previous work these collisions will be hard 
or impossible to observe, and we will exclude them from 
our analysis, while noting that they still comprise some- 
thing of an open problem as to precisely how this infinity 
should be regulated. 

A second sizable contribution to V^ obs results for large 
Hp I Hi from bubbles appearing on the conformal dia- 
gram near future infinity of the false vacuum. This en- 
ters the observer's past lightcone when the "hat" in the 
conformal diagram (Fig. [TJ extends high enough, which 
in turn happens for Hp 3> Hj. This effect is largely in- 
dependent of £0 and corresponds to "late time" collisions 
of all angular scales, and nearly uniformly distributed in 
cos 9. 

Note that many late-time collisions will result in effects 
on scales greater than the observer's full sky. To see why 
this is the case, we can coordinatize the exterior dS in 
terms of angular scale and z c . This was done in Ref. [23] 
which showed that in the observation frame, each bub- 
ble intersecting part of the observer's LSS (but not en- 
compassing the entire surface) with solid angular scale 
ip and dS-invariant distance corresponding to z c can be 
mapped uniquely to a nucleation point (r] n , T n ) in the ex- 
terior dS. A straightforward generalization of Appendix 
D of Ref. [33J allows us to exchange the angular coordi- 
nate ip for one describing the size of greater than full-sky 
bubbles as well; this quantity, which we call p, is defined 
via 



■1 



sinh-ff/Tis sinh 6. 



2 „2\ 



2H lZcl 



[1-H]z 2 c ) cashHms 



+ v sinh HjTig cosh £;s 



(7) 

where 6s = J dr /a(r) is the comoving radius out to 
which the observer can see on the LSSQ For -1 < p < 1, 



p corresponds exactly to —cosip/2. The coordinate p 
effectively measures the distance out on the observer's 
last scattering surface to its intersection with the future 
light cone of a given collision in the observation frame. 
Values p < — 1 correspond to late-time collisions which 
do not intersect any of the observer's LSS, while p > 1 
corresponds to a greater-than-full-sky collision. 



The foliatiorf 
shown in Fig. p 



of the exterior dS in terms of (p, z c ) is 
_for H F I Hi = fO, 6s = 1 on the left, 
and for Hp/ 'Hi = 100, 6s = 0.05 on the right. The col- 
ored solid lines correspond to p values such that \p\ < 1, 
i.e. nucleation centers of bubbles with angular scales 
< ip < 2tt. The (0-dependent) initial value surface 
is not shown. Notice that for small values of 6s (as re- 
quired by observation) and Hp /Hi » 1, bubbles of vir- 
tually all angular scales (including p > 1) receive a large 
4- volume contribution from points near future infinity of 
the exterior dS space. This is true regardless of the po- 
sition of the initial value surface and hence of the angles 
9, <f>. Thus, neglecting early-time bubbles and for values 
of Hp I Hi, 6s we are interested in, an observer should 
expect to see a nearly scale-invariant, isotropic distribu- 
tion of collision bubbles, regardless of his or her position 
in the observation bubble. (This typically includes many 
greater-than-full-sky late-time bubbles, as is evident from 
Fig.[2j) As the relevant constant- p curves are compressed 
into a very thin region in this case, from Fig.[2]one would 
expect the scale-invariance to hold true independent of 
z c , although the distribution of bubbles itself will have 
z c -dependence. 

Explicitly calculating the distribution of bubbles over 
angular size by integrating the 4- volume element over the 
constant-p contours in the relevant region yields 



dN 
dpdfl 2 



XHZ 



El 

Hi 



6. 



(8) 



for 6s - > 0, which is indeed flat in p. This distribution 
was first obtained by Ref. [IT] and dubbed the "disinte- 
gration of the persistence of memory" due to its isotropy. 

Integrating the 4-volume element over all values of 
p > — 1 and z c in the region depicted in Fig. [l] yields 
the expected total number of bubbles intersecting the 
observer's LSS, given by 



N: 



4ttA 

zhJ 



6, 



(9) 



for observers at 6) ~ ^ 00 Q Limiting the integration to 
— 1 < p < 1 removes greater-than-full-sky bubbles from 



4 Measurements of the curvature density [231 125] dictate that 



2\/n^ < 0.18. In what follows, we thus assume £; s is small, 
typically taking £; s i=s 0.05 as an illustrative value. 
For p > 1 there can be two points in the false vacuum cor- 
responding to the same value of (z c ,p) as seen on the LHS of 
Fig. [2] However, we are most often interested in \p\ < 1, for 
which the foliation is one-to-one. 

Virtually all observers, except a set of measure zero, are expected 
to fall into this category. 
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FIG. 2. Foliation of a 9 = conformal slice of the external de Sitter spacetime for an observer at £o = with lines of constant 
p, z c for the (unrealistic) values £; s = 1 and Hf/Hi = 10 on the left, and for £; s = .05 and Hf/Hi = 100 on the right. The 
observer's LSS is sketched by the dashed black line in the observation bubble (this curve will be compressed up towards the 
tip of the observation bubble for the case shown on the right). Solid black curves correspond to different values of z c , with the 
vertical line in the center corresponding to z c — Hp 1 in both cases. The dotted curves represent different values of p, with 
p = — 1, 0, 1 denoted by the solid green, red, and blue curves, respectively. The effect of decreasing £; s and increasing Hf/Hi is 
shown on the right. The same range of p values are plotted, only now all constant-p curves depicted on the left are compressed 
into the thin colored region. 



the count and yields 



H„ and falls off as 



for large 



N ohs 



16ttA (El 



3Hp 



Hj 



'tin 



(10) 



Comparing N to iV D b s we see explicitly that for each 

— 1/2 

late time sub-horizon bubble, there will be ~ £o^e 3> 
1 full-sky bubbles to account for when computing the 
effects of bubble collisions. As we will argue in the fol- 
lowing sections, however, collisions that encompass the 
entire visible portion of the LSS will result in adiabatic 
superhorizon perturbations to the gravitational potential 
at the LSS and hence not contribute significantly to the 
CMB temperature spectrum. 

From Eq. [lOj there are many possible potentials and 
nucleation rates giving rise to any particular value of 
^obs- Our approach will be to treat N \, s as a free pa- 
rameter when computing the expected CMB tempera- 
ture spectrum and we will not attempt to address how 
feasible or likely scenarios with a given iV bs might be, 
asi de f rom the following general comments. First, from 



Eq. 



10 



even for \H F A 



<C 1, one can still expect 7V obs ^> 1 



as long as Hp /Hi is large enough to compensate for the 



small nucleation rate. Second, it is clear from Eq. 10 
that iVobs — 1 is likely a very special case, requiring the 
precise balance between features of the potential through 
Hp, Hi, f2 c , and A as mentioned previously. 

Finally, one can also compute the distribution of col- 
liding bubbles with respect to z c . This calculation was 
carried out by the authors of Refs. pj [53] who found that 
the distribution dN/dz c dQ2 is peaked sharply around 



III. THE PROFILE: SINGLE BUBBLE CASE 

Calculating the CMB temperature spectrum resulting 
from a multiple collision scenario requires an understand- 
ing of the temperature profile associated with a single 
bubble collision, which we discuss here, as well as how 
to combine many such collisions, which is discussed in 
the next section. For more detailed derivations of the 
single-collision temperature profile we refer the reader to 

Refs. pa era eg. 



The Single-Bubble Profile As A Small Potential 
Perturbation 



As noted above and in previous studies, the symmetry 
of a bubble collision leads to an azimuthally symmetric 
perturbation on an observer's LSS, filling a disk that con- 
stitutes the triple-intersection of the last-scattering sur- 
face, the observer's past lightcone, and the future light 
cone of the bubble collision. Within this disk, given the 
necessary slow-roll inflation in the observation bubble 
and the lack of large anomalies on the actual CMB, the 
effect of the collision is necessarily small, and can be de- 
scribed as a perturbation to the gravitational potential 
&(x) at the reheating time T r /j. In conformal Newtonian 
gauge and retaining only the term linear in slow roll (sec 



G 



Ref. [T7]), the gravitational potential resulting from the 
inflaton perturbation is given by 

lb lb 
${T r h,x) = a^ r h(cos6 — cos — )0(cos# — cos — ) (11) 

where coordinates are chosen so that the collision is cen- 
tered around 9, (j> = with 9 = ip/2 the causal boundary 
of the collision, Q(x) is the Heaviside step function, and 
£ r h is the comoving distance at which the observer's past 
light cone intersects the reheating surface. The param- 
eter a contains all the information about the inherent 
strength of the perturbation, which depends on the kine- 
matics and the shape of the inflationary potential. 



Equivalently, the gravitational potential Eq. 11 can be 
exchanged for the gauge-invariant curvature perturba- 
tion, which at late times during inflation is simply £ oc $; 
the proportionality constant here depends on the infla- 
tionary potential and a. Once £ is known, one can use the 
Sachs- Wolfe approximation [26] AT/T oc ((x,ti 8 ) to de- 
termine the CMB temperature spectrum resulting from 
the collision. Here the factors relating £ to AT/T de- 
pend only on the background cosmology. Absorbing a 
and the other relevant parameters fixed by the poten- 
tial and intra-bubble cosmology into the parameter a, 
and neglecting (temporarily) the inflationary tempera- 
ture fluctuations, as well as the evolution of the pertur- 
bation between reheating and decoupling^ we arrive at 
the approximate CMB temperature profile for a single 
collision 



T(6, 4>) = T [1 + a£i„(fj, - Po)Q{v - Mo)] 



(12) 



where p = cos 9, po = cos ijj/ 2 for a collision centered at 
the north pole and To is the average temperature of the 
unaffected portion of the LSS. 

In sum, the form of the perturbation in the current ap- 
proximation depends only on the residual SO (2, 1) sym- 
metry of the post-collision spacetime, and all the infor- 
mation about the microphysics of the collision, including 
the details of the underlying potential and background 
cosmology, is encoded into a single parameter a describ- 
ing the inherent "brightness" of the collision on the CMB. 
This fact will simplify the calculation of the expected 
CMB temperature spectrum for a multiple collision sce- 
nario in Sec. IVl 



B. The Frame Shift 

The derivation of the temperature profile Eq. [12] is 
valid for scenarios in which the observer is born comoving 



7 Technically one should take into account the evolution of Q be- 
tween reheating and decoupling by using the appropriate transfer 
functions. This was carried out in Ref. |17l which found that the 
results are in good agreement with Eq. |12| (This is expected, 
since for Hp /Hi S> 1 the conformal time elapsed between re- 
heating and decoupling is negligible.) For the remainder of this 
study we thus approximate the reheating time T r h ~ r; s . 



with respect to the unperturbed portion of the last scat- 
tering surface. In this case, the perturbation from the 
collision enters the past lightcone of the observer at late 
times. However, it was shown in Ref. [5] that 0(3, 1)- 
invariance can be spontaneously re-generated in the re- 
gion to the future of a collision, thus "just as many" ob- 
servers can potentially form comoving with the perturbed 
portion of the last scattering surface. (Such observers 
were dubbed "foreign-born" in [53] .) As described be- 
low, these observers are typically those that would see a 
bubble with p > (covering more than half of the sky) . 
By the arguments of the last section, most observers who 
witness any bubble collisions at all will see many with 
p > 1, and hence be of this type. 

We argue, however, that the effects of p > 1 bubbles 
will be largely invisible, and that the effects of p > 
bubbles will be degenerate with (and convertible into) 
those of p < bubbles. To see this, we can adapt the 
arguments set forth in Ref. [37] as follows: in a ACDM 
universe, the total temperature anisotropy in a direction 
n on the sky is given by a sum of the Sachs- Wolfe (in- 
cluding integrated effects) and Doppler anisotropies in 
conformal Newtonian gauge 



AT 

AT 

~T~ 



(A) 
(n) 



- $(r ;s ,n^ s ) - 

SW+ISW 

v(t ,0) - v(ri s ,h£ ls ) 



-t-eV) 



(13) 
(14) 



where the velocity v induced by a superhorizon pertur- 
bation is given by 



v(t, x) ~ V<f> 



d\na 



(15) 



(note that the we have neglected the curvature of the last 
scattering surface, a good approximation for £; s -C 1). 
Since the above expressions for the temperature differ- 
ence are linear in $ and its derivatives, to lowest order in 
velocity and considering only the linear piece of the 
various contributions for a gravitational potential of the 
form $ = + $ 2 are simply the sum of the individual 
contributions: 



AT AT. 



AT 



(16) 



Consider a particular collision intersecting the ob- 
server's LSS with angular scale tp- Neglecting the cur- 
vature of the last scattering surface, one can imagine the 
collision inducing a planar perturbation 



= a(x ~ X int )Q{x - X int ) 



(17) 



at T; s , where both nucleation centers lie along the £-axis, 
and Xi n t denotes the intersection of the future light cone 
of the collision with the LSS. One can rewrite Eq. 17 as 
a superposition of two modes &i(x), &2{x) as 

= a(x-x int )+a(x int -x)Q(x int -x) = ^ l (x)+^ 3 (x). 

(18) 



Here <5>i(x) is a superhorizon mode while $2(3?) is equiv- 
alent to that of a collision centered around the antipodal 
point of the original collision on the celestial sphere with 
angular scale ifi' = 2ir — ip and strength a. 

A foreign born observer emerging from the last scat- 
tering surface at x > Xi n t that is initially comoving with 
the perturbed LSS only sees the superhorizon mode $1 
at early times, with $2 only entering his or her past light 
cone at late times. As a result, the contribution of $1 to 
the observed CMB temperature spectrum will vanish as 
in Ref. [37] and using Eq. 16 the resulting CMB temper- 
ature fluctuation will be 
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AT 



2 ■ 



(19) 



The effect of a collision on the CMB temperature profile 
for a foreign-born observer from a collision centered at 
(^Oj'/'o) with angular scale ip is therefore equivalent to 
the profile seen by a native-born observer for a collision 
centered at the antipodal point (it — 9q,2tt — 0o) with 
scale ip' = 2ir — ip; the observed perturbation is shifted 
by going between the native and foreign born observers' 
reference frames. 

The above argument applies to observers whose past 
worldlines intersect the perturbed portion of the last scat- 
tering surface for a given collision. This will clearly be 
the case for collisions encompassing the observer's entire 
LSS (p > 1), whose perturbations consequently vanish 
in our approximation. This allows us to neglect the po- 
tentially large number of full-sky collisions accompany- 
ing A^bs > 1 in calculating the CMB temperature power 
spectrumj^] 

This "frame shift" will also tend to apply to collisions 
with p < 1 that cover more than half of the observer's 
sky. For a native born, comoving observer originating 
from x ~ Xi n t, the collision will nearly bisect the LSS. 
This situation is depicted by observer A in Fig. [3j As 
long as a foreign-born observer's worldline does not inter- 
sect those of native comoving x ss Xi n t observers (which 
should typically not occur, since late time peculiar veloci- 
ties will be dominated by the effects of processes uncorre- 
cted with the collision itself), the foreign born observer 
will see a greater-than-half sky collision (observer B in 
Fig. [3]) . Conversely, native-born observers will only see a 
collision with p > if they also acquire a significant pecu- 
liar velocity in the direction of the collision (c.f. observer 
C in Fig. [3]). The farther away from Xi n t the observer is 
born, the larger the peculiar velocity will need to be for 
him or her to observe a greater than half-sky bubble that 
is not frame shifted. Thus we expect that most p > col- 
lisions will result in frame-shifted profiles, appearing to 
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FIG. 3. The universe with one bubble collision at decoupling. 
We show three different observers who see a collision with pro- 
file Eq. 
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taking up (A) half of the sky, (B) more than half of 
the sky, and (C) less than half of the sky. Their present-day 
last scattering surfaces are shown by the blue circles, corre- 
sponding to the intersection of their past lightcones with the 
decoupling time slice. Unless the observers have acquired a 
large peculiar velocity along the x direction, the last scatter- 
ing surfaces at some earlier time are approximately given by 
the dashed blue lines for each observer. In this case, observer 
B is foreign born and sees a frame shifted collision while C 
does not. Observer A is a boundary case and could have 
emerged from either portion of the LSS. In our approxima- 
tion, we take all observers born with x > Xi n t to observe a 
frame shifted collision 



take up less than half of the observer's sky, while most 
with p < will be unaffected. 

For illustrative purposes we will work under the as- 
sumption that all greater-than-half sky bubbles corre- 
spond to frame-shifted collisions. We implement this in 
the bubble profile by mapping all such collisions to less- 
than-half sky events centered around the antipodal point 
of the original nuclcation center on the sky. Relaxing this 
assumption does not significantly affect our conclusions, 
since this prescription primarily effects the dipole, reduc- 
ing it compared to the case where the frame shift is not 
taken into account. 



IV. MULTIPLE COLLISIONS: GENERALIZING 

TO iVobs > 1 

A. A Multiple Collision Scenario As The Sum Of 
Its Parts 



It is possible that there are other observable effects associated 
with collisions on these scales. We do not consider this possibility 
here and in the remainder of this work simply neglect bubbles 
with scales p > 1, cautioning the reader that the compatibility 
of some large- Afobs scenarios with observation may be spoiled by 
the inclusion of such effects. 



We wish to generalize the results in Sees. [TTJjTTT] to the 
case of A Q bs > 1- Colliding bubbles that do not overlap 
in their influence prior to the last-scattering surface can 
simply be treated as independent. Events where bub- 
bles do overlap prior to this time can divided into "pre- 
collisions" occurring before either bubble intersects the 
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observation bubble, and overlaps occurring within the 
observation bubble, prior to last-scattering. We will first 
show that the first category are uncommon and hence 
generally unimportant, then discuss the second category. 

Consider the disturbance from a pre-collision impact- 
ing the observation bubble wall at To in the observation 
frame. We can model this disturbance as a shell of radi- 
ation emanating from the pre-collision event. Observing 
this radiation shell implies that one of the iV bs collid- 
ing bubbles, which would have impacted the observation 
bubble at T , was intercepted by an intervening bubble 
before it could do so. For a given bubble with nucleation 
centeiP](z c , p) colliding with the observation bubble wall 
at To , we can ask how many intervening bubbles of a par- 
ticular type it expects to encounter along the way to the 
observer's bubble. We denote this quantity as N int (z c , p); 
for a fixed A, it is given by 

N iat (z c ,p) = XX(z c ,p) (20) 

where a conformal slice of the 4-volume T{z c ,p) is de- 
picted by the region shaded yellow in Fig. [4] Note that 
iVi n t includes bubbles with p > 1 in its count, and that 
these bubble need not nucleate at the same (9, 4>) values 
as the collision bubble. 

The 4-volumc T{z Cl p) corresponds to the portion of 
the past light cone from the would-bc intersection point 
of the collision bubble with the observation bubble wall, 
(jio,T ), not also within the past light cones of (z ci p) 
and the origin. Intervening bubbles nucleated outside of 
the past light cone from (770, To) cannot produce radia- 
tion from the pre-collision observed at To, while requiring 
that the collision and observation bubbles were not born 
within the intervening bubble eliminates the 4-volume in 
the past light cones from (z c ,p) and the origin, respec- 
tively. Defined this way, I(z c , p) will not include pre- 
collisions whose radiation enters the observation bubble 
farther up the H2 collision surface (i.e. T > To), but 
these events will be counted in iV; nt for the correspond- 
ing intervening bubble. Thus, for a given -/V b s , we can 
look at iVi n t for each observable bubble and if this quan- 
tity is small in all cases, we can safely neglect the effects 
of all pre-collisions on the effective temperature profile. 

Using the 4-volume element in Eq. [3j we can estimate 
iVi n t for a given collision bubble nucleated at (z c , p). This 
calculation is detailed in Appendix [XJ The results are 
illustrated in Fig. [5] where A is determined by requiring 
^obs = 1- For large Hp /Hi, the curves of constant-/? for 
a given z c are compressed as in Fig.[2]and so the curves in 
Fig. [5] are nearly identical for all p in the interval [—1,1]. 

From Fig. [5] we see that most bubbles will indeed not 
collide with an intervening bubble. This is because al- 
though many bubbles can nucleate as T — > tt/2, very 
few of them will have a chance to intersect the collision 



9 As discussed in Appendix [X] we are free to choose 8 C = <f> c = 
without loss of generality. 



bubble wall before it impacts the observer's bubble. The 
past light cone from the intersection point of the collision 
bubble with the observation bubble wall contains very lit- 
tle 4-volume near the future infinity of dS, only opening 
up to enclose more 4-volume for smaller T. As a result, 
regions for which dV± starts to diverge do not contribute 
very significantly to I. Note also that smaller values of 
z c correspond to smaller values of N- mt . This is because 
for large Hp /Hi and the small £; s we consider, the colli- 
sion bubble nucleation points corresponding to small z c 
are very close to the observation bubble wall, resulting in 
a smaller 4-volume T(z c ,p) and hence smaller A^ nt . For 
larger z c , this cancellation is not as severe and -/Vj nt can be 
close to 1 for small Hp /Hi, however this will necessitate 
a large A which is in tension with the requirement for eter- 
nal inflation. Also, increasing iV bs will shift the curves 
upward, however for cases of interest with Hp /Hi 1, 
£; s <C 1, we will still expect N int <C 1 unless iVobs is very 
large as well. Thus, in what follows we can safely neglect 
any interactions between the bubbles before they impact 
the observation bubble wall. 

Although given our parameters of interest, bubbles 
only rarely "pre-collide" , they should commonly overlap 
within the observation bubble if A^bs 3> 1, especially 
considering p > 1 ("all sky") collisions. If these over- 
laps occur after the effects of the collisions have entered 
the perturbative regime, the various perturbations will 
simply superpose. This will not necessarily be the case 
for large iV bs, however. Consider the 2-sphere defined 
by the intersection of the observer's past light cone with 
the observation bubble wall. Given a uniform distribu- 
tion of collision bubble nucleation centers, the points of 
intersection of the collision bubbles with the S2 will be 
distributed with average distance ~ ttRq /iV b s separating 
any two of these points, where i?o is the observation bub- 
ble radius at the time of intersection. Since Rq ~ HJ 1 , 
one would expect the future light cones from these points 
to intersect early on for large N \, s , and so in principle 
the effects of the various collisions may interact with each 
other before enough inflation has occurred to dilute the 
disturbances to the point of being perturbative. 

To study these interactions would require details of the 
microphysics governing the collisions, as in the numerical 
study of Ref. [21] . Such considerations are beyond the 
scope of this paper, and so in what follows we neglect 
these interactions: in treating the effects of the collisions 
inside the observation bubble, we approximate the vari- 
ous disturbances as simply superposing on one another. 
The resulting effects on the CMB temperature will then 
be described by the sum of Sachs- Wolfe contributions 
from the individual collisions (including the frame shift). 
We intend to test this approximation in future numeri- 
cal studies, but for the time being we content ourselves 
with this treatment, as we expect it to characterize (at 
least approximately) the effects of overlapping collisions 
on the CMB. 

With the above assumptions, a multiple collision sce- 
nario can indeed be analyzed as a collection of single col- 
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FIG. 4. A conformal slice of the 4- volume available to nucleate an intervening bubble for a collision bubble nucleated at (z c , p) 
(with 9 = (j) = 0). The region I(z c , p) is given by the 4-volume within the past lightcone from the point on the observation 
bubble wall with T = To (in green), excluding the the 4-volume within the past lightcones from (z c , p) and the origin shown in 
red. 



iVtat 




H F 



10 100 1000 10 4 10 5 10 6 Hi 

FIG. 5. The expected number of intervening bubbles as a 
function of Hf/Hi for one collision bubble (iV b s = 1) nu- 
cleating from z c — 0.5,1,5 (green, blue, red) and p = 1 in 
the false vacuum and with £; s = .05. Varying p in the in- 
terval [—1,1] has little effect on the curves, as expected from 
the discussion surrounding Fig. [2] Increasing JV b s will shift 
these curves upward, however for Hf/Hi > 1 we still expect 
JVtai « 1- 



be written as 

T(n) =1* (1 + A(n)) 

X [l + ^ a i^s(/"i ™ COS y )8(/l; - COS y) 

+ a iiis{-^i + cos + COS y ) 

N 2 

( 21 ) 

where we have defined the angular distance for a given 
collision as 

fii = sin 8 sin 9^ (cos cj) cos 4>oi + sin <p sin 4>m ) +cos 9 cos #oi 

and where (^oii'/'oi) is the angular position of the nu- 
cleation center of the ith bubble on the observer's sky, 
Ni (N2) is the number of bubbles affecting less than 
(more than) half the sky, is the inherent brightness of 
a given collision, Tq is the average temperature neglecting 
any perturbations, and we have included the inflationary 
density perturbations A(n) which we approximate to be 
unaffected by the collisions. Note that in order to treat 
effects of the collisions on the CMB temperature pertur- 
batively as we have done, we should have -/V b s ai£is <C 1 
for all collisions. With a few additional assumptions we 
can use Eq. [21] to extract some general features of the 
effects of multiple collision scenarios on the CMB. 



B. Additional Assumptions 

lision events. Accounting for the frame shift (and thus 

neglecting bubbles with p > 1), the resulting CMB tern- Several more useful observations and simplifications 

perature profile for a given multiple collision scenario can can applied to the picture above: 
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• We neglect the possibility of (p < 1) bubbles form- 
ing within other (late time) bubbles. For small A 
this is reasonable, for similar reasons as the neglect 
of pre-collisions: the overlap between the past of an 
observer's last-scattering surface and the future of a 
given bubble collision is generally of order Hp , so 
the number of nucleations from that region, which 
would consitute bubbles-in-bubbles, is small unless 
A > Hp*. Arguments along these lines have been 
realized in the past (see e.g. Refs. [6] fTTj) . 

• We assume that only one decay channel is rele- 
vant for nucleating colliding bubbles from the par- 
ent vacuum. For a long-lived vacuum this does not 
seem unreasonable: since A = Ae~ SB , where Se is 
the Euclidean action of the CDL instanton and A 
is a prefactor containing quantum corrections, field 
trajectories resulting in small differences in Se will 
tend to engender rather large differences in A so 
that one dominates. For a recent discussion along 
this line of thought, see Ref. [28]. This assump- 
tion allows us to consider only one type of collision 
bubble relevant for the multiple collision scenario. 

• We imagine that the potential V(<p) in the neigh- 
borhood of <p after any of the collisions is relatively 
featureless, i.e. V(<p + S<pi) « V(p + Stpj) and like- 
wise for V . Since we require a period of slow roll 
inflation to the future of all collisions, which implies 
a flat potential near the post collision field values, 
this assumption is typically satisfied a posteriori. 
We note that it is difficult to see how this assump- 
tion, or even that of a multiple collision scenario, 
could be satisfied for small-field inflationary mod- 
els (i.e. models where the width of the inflationary 
trajectory Ap < M p i) without very small values 
of Sip, since even small perturbations will tend to 
cause the inflaton to overshoot the inflationary re- 
gion of V(cp) [21]. 

The main consequence of the above assumptions for 
our purposes is that the brightness parameters (Xj only 
depend on z c and p, by the following reasoning. For a 
fixed potential with colliding bubbles of only one type, 
the strength on of the perturbation at r rn depends on 
the kinematics through Sip. In the collision frame, Sip 
is determined only by the parameter z c , while the boost 
back into the observation frame generally depends on p. 
However, since we are most interested in cases where 
He I Hi 1, collision bubbles with p in the interval 
[—1, 1] arc all mapped to a very small range of r\, T for a 
given z c , and hence all receive virtually the same boost. 
Also, a is related to Sp by factors which depend on the 
background cosmology as well as V(ip), V'(<p), but by the 
third assumption above, this factor should be similar for 
all collisions. Consequently, a ~ ot(z c ) and the brightness 
parameters in Eq. [21] in our approximation depend only 
on the kinematic variable z c and the underlying poten- 
tial, up to some overall normalization. We will exploit 



this fact in taking observer expectation values for the 
temperature spectrum in the following section. 



V. EFFECTS ON CMB TEMPERATURE 
POWER SPECTRUM 

Given that, in many cases of interest, the effects of 
multiple bubble collisions may be approximated by su- 
perposition, it is worth considering the observational con- 
sequences of such scenarios by generalizing previous work 
to the case of many bubbles. Armed with the tempera- 



ture profile Eq. 21 the remainder of the paper constitutes 
a first look at multiple collision signatures, by computing 
their effects on the CMB temperature power spectrum. 

Our setup and framework follows that suggested in 
Refs. pUl [29], with the appropriate generalizations to 
multiple collisions. Denoting the sky-averaged observed 
temperature as To, we expand the observed temperature 
fluctuations for a given set of collisions in terms of spher- 
ical harmonics 



ST(n) 



T-T 



= ^Tt lm Y lm {h). (23) 

l, in 



Wc define 



A{h) = f{l + w{n))-l (24) 
B(n) = A(n) (25) 



JV, 



+ ^ 0!i^i s (-pi + cos y )6(-M* + cos y ) 

iVj 

(26) 

where / = Tq/Tq, and A(n) is the standard inflationary 
temperature fluctuation in direction n. (We assume To = 
2.725 K throughout our calculations.) 

Using these definitions we can recast Eq. [23] as 



ST(n) 

Tn 



= A(n) + f(l + w(n))B(n). (27) 



Like the temperature contrast of Eq. |23[ we can expand 
the functions w, A, B in spherical harmonics with cor- 
responding coefficients Wi m , ai m , bi m . Then, using the 
properties of products of spherical harmonics, the coeffi- 
cients ti m are given by 

U m = a lm + fb lm + fJ2 J2 w hmA 2 m 2 R l C ll2m2 ( 28 ) 

where the R terms are Gaunt coefficients, given in terms 
of Wigner 3 — j symbols by 



R i imihm2 =( 1)m / (2i 1 + l)(2/2 + l)(2i + iy 



'1 '2 ' \ / n I2 I 



(29) 



/ \ mi rri2 —m 



11 



The conventions we use here are those found in Ref . [3D] j 
to which we refer the reader for additional details. Note 
that the isotropy in the corresponding expressions in 
Ref. [10] is broken as a result of including more than 
one collision. 



A. The Two-Point Function 

We are interested in the expected angular temperature 
power spectrum for a CMB sky affected by many bub- 
bles. To compute this quantity, one must evaluate the 

using Eq. 



two-point function (ti m t* m ), using Eq. 28 to average over 
both the ensemble of density fluctuations from inflation 
(hereafter "ensemble averages"), as well as over the dis- 
tributions governing the various collision parameters. 



1. Ensemble Averages 

First let us perform the ensemble averages. Consider 
a given set of collisions B arising from a fixed potential 
and -/Vobs- Following previous work, we approximate the 
inflationary density perturbations as being unaffected by 
the collisions. The only terms affected by the ensemble 
average for a given set of collisions are those containing 
bi m in Eq. |28| Since the density perturbations are as- 
sumed to be Gaussian with zero mean, the one point 
functions (bi m ) vanish for any given set of collisions 
and only terms with two-point functions (bi m bf m ) will 
survive (here the brackets () denote ensemble averages 
with B fixed) . If the density perturbations are unaffected 
by the collisions, the corresponding two-point function is 
given by 



(hmb* r , 



C bb 



(30) 



where C\ h is the correlation function in the absence of 
any collisions (note that (bi imi b^ 27n2 ) B — for l\ ^ 
We obtain the C\ h from CMBEASY [3T] using concordance 
WMAP7 values [32] for the relevant cosmological param- 
eters. Performing the ensemble average for this set of 
collisions yields 

(ti m t*i m ) B =ai m a* lm + f 2 (c bb + 

J2w hmi R l l ^ lm Ct b + cc.+ 



wi imi w l3m3 n lm n lr 



C, 



(31) 

where the sum in the third line is over i = 1,2, 3. From 



Eq. 31 we can calculate the CMB angular temperature 
power spectrum Cf for a set of collisions B via 



C? = 



21 



(32) 



2. Averages Over Collision Scenarios 

To understand the generic features of a multiple colli- 
sion scenario, we wish to average Eqs. [31] [32] over the set 
of possible collision scenarios, {B(N a \ is )} with the under- 
lying potential and N ^ s fixed. This yields 



(tlmtfm) = (ai m a* lm ) + (f) (c bb + 
{wi imi )H lm L l 



cc. 



El... „,,* \ nl 3 m 3 l 2 ra 2 r>l 1 m 1 l 2 m 2 fibb 



(33) 

where the brackets now denote averages over {B(N b s )} 
for a given potential and N Q ^ S . From this, one can com- 
pute the expected CMB angular temperature power spec- 
trum via 

(Cl) = Mm) ■ (34) 



B. Performing The Averages 

To perform the averages over the various B(N \ 1S ) in 
Eqs. 33| 34 requires averaging the one- and two-point 
functions of Wi m over the distributions of collision pa- 
rameters, which we now address. 

Consider an observer with N Q ^ S collisions intersecting 
the visible portion of the LSS. Any one of the -/V b s col- 
lision bubbles can be thought of as being selected out of 
a probability distribution V(z c , if), 9, 4>) that corresponds 
to the distribution over nucleation events (as discussed 
above, we will also assume that this suffices to determine 
the overall effect of all collisions). Neglecting the effect 
of bubbles in bubbles (by effectively allowing bubbles in 
bubbles that superpose) renders each nucleation event 
statistically independent and so V should be the same 
for each bubble, given by normalizing dN / dz c dpdD.2 to 
1. For large Hp /Hi and small the 4- volume element 
dVji{z c , p, 9, (f) at a given z c is very nearly independent of 
p, 9 and </>. Thus we can write 



dN(z c , p, i 



XH F 4 j(z c )dz c dpdQ 2 , 



(35) 



where j(z c ) is a known function of z c which one can com- 
pute from the definition of z c . From the discussion in 
Sec. |nj j(z c ) is peaked around z c ~ Hp 1 , and decreases 
for larger values of z c . Since we are concerned only with 
bubbles such that \p\ < 1, we change variables back to tjj 



and normalize Eq. 35 to yield 



1 ip 
V{z c ,il),9,4>) = — j(z c ) sin-. 

167r 2 

Here j{z c ) denotes the normalized j(z c ). 



(36) 
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FIG. 6. Left: The dimensionless quantity Wi which enters into our final expression for (Ci), plotted for < I < 40. The I = 0, 1 
values are reduced relative to the case not accounting for the frame shift, wi encodes the shape of the spectrum resulting for 
the isotropic and scale-invariant dist ribution of p; its contribution to (Ci) is multiplied 
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plotted for 2 < I < 40. Ri also enters into Eq. 
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denote the quantity Ri from Eq 

that the dips for I = 2, 40 are a result of excluding h < 2 and h > 40 in the sum in Eq 
for a concordance ACDM cosmology from CMBEASY are shown by the large black points. 



by ty N obs a2^f s . Right: The red points 
carrying a factor of n 2 N b B ct2£,i s - Note 
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The theoretical values Cf b obtained 



The probability density V governs a single collision, 
but we can use it to average over particular realizations 
of multiple collision scenarios with -/V b s fixed. This is be- 
cause the various quantities entering into Eq.[33]consist of 
the sums of contributions from 7V b s individual collisions 
and, if each nucleation event is statistically independent, 
the averages over {B(N D h s )} and the sums over bubbles 
commute. Thus, the expectation values of the one- and 
two-point functions for many bubbles are the sums of the 
corresponding expectation values for individual bubbles, 
which are averaged over V . For example, the expectation 
value of the one-point function wi m for a scenario with 
-Nobs collisions is given by 



(Wim) = ( W lm 



(wj m )=N ohs (w? m ) (37) 



where (wf^ is the expectation value for a single collision, 
which is the same for all collisions because of statistical 
independence. Similarly for the two point function, 



( 



w h mi w hm2 



> = (£' 



'l imi w l 2 m 2 



E 



*0 \ 

2"J2 / 



+ N ohs (N ohs - 1) « roi ) (w*° 



(2TO2 



> 



We can then use the distribution V to perform the single- 
bubble averages in the last equalities of Eqs. [37j |38j The 
details of these calculations can be found in Appendix [5] 
Consider first the averages over z c . As discussed in 



Sec. IV B the dependence of the temperature profile on 
z c for Hp /Hj ^> 1, £i s <C 1 arises only through a, 
which additionally depends on the underlying potential 
and microphysics but is independent of ijj in this limit. 



Therefore we can absorb all of the uncertainty and model- 
specific information pertaining to the averages of a over 
z c into two parameters, 



oti = (a) = J a(z c )j(z c )dz c (39) 
a 2 = (a 2 ) = / a 2 {z c )j{z c )dz c . (40) 



Doing so allows us to treat 0^2 (along with N b s ) as free 
parameters and work out (C;) in terms of these quan- 
tities; we defer a study of realistic ranges of ai.2 to fu- 
ture work. As a consequence, our computed (Ci) will 
describe the expected power spectrum in terms of iV bs, 
Q!i, 2) assuming that a scenario described by these values 
can occur. Recent numerical studies [H] suggest that the 
strength of the inflaton perturbation from a collision (and 
hence a) are primarily dependent on the underlying po- 
tential, with only a mild dependence on the kinematics. 
If this is the case, then for a fixed potential the averaged 
parameters a\, 0*2 should provide a good characteriza- 
tion of the expected signal for most observers. Note that 
by Jensen's inequality, a.2 > a 2 . 

Performing the averages over V , the expectation value 
of the one-point function (wi rn ) is given by 



(Wl m ) = — iVobsai&a<5j0<W>- 

6 



(41) 



Only the monopole contributes to the one-point functions 
by the orthogonality of the spherical harmonics. This 
monopole piece enters at all Z-values by multiplying the 
contribution of the inflationary perturbations in Eq. |33| 
Moving on to the two-point function, and motivated 
by the discussion in Appendix IbJ we define the quantity 
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wi as 



cring the last term in Eq. 33 we define 



Wl 



2 & 

\wi'0'(ip)\ sin— 5u>dip 



where 



Y, 



I'O' 



1> 



(42) 



dx (tt — 



YJq, (cos 1 x) ( x + cos — J dx9 — 7r) 

)B $ V 2 / 

(43) 

As detailed in Appendix [B] primed indices as in Eq. [43] 
signify quantities evaluated in the polar frame, in which 
a collision is centered at the north pole of the celes- 
tial sphere. The quantity wi isolates the /-dependence 
of (wi imi 'Wi 2m2 ). We plot w t for multipoles I < 40 on 
the left in Fig. [6j In terms of wi and neglecting the 
monopole contribution to a/ m , the two point function 
(ai imi a*i 2m ,) = (f 2 ) (m imi Wi 2m2 ) where ( see Appendix 



Km, < m2 ) 



— 2 j - 1 — w h d h i 2 d mim2 + 



A bs(Aobs — 1) ^ <^iO < ^20 < Wo^m 2 

The first piece contains the expectation value of the in- 
dividual two-point functions and the second piece cor- 
responds to products of the one-point functions. This 
contribution is statistically isotropic. 

I t re mains to average the Gaunt coefficient piece of 
Eq. 33 as well as to compute the quantity (/ 2 ) . Consid- 



m Zi,mi,/2,^: 



2Ji + l 



R 



l 1 m 1 l 2 m2 
Im 



C\ 



1,1, 



(45) 



We plot Ri on the right in Fig. along with C\ h for 
reference. Clearly, R\ <C w; since by Eq. [45] it involves 
the product of perturbation terms w;C ; h& , and the Gaunt 
coefficients are typically O(10 _1 ) or smaller in our com- 
putations. In addition to i?z, there is also a contribution 
^ 2 afo 2 s i?°(iV obs ) to d where 



^(AW) = 



jVobs(jVobs - 1) 
1447T 2 



C, 



bl, 



(46) 
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This term arises from the second line of Eq 
inserted into Eq. [33[ 

To compute (/ 2 ) we must take the ensemble average of 
the inverse of the sky-averaged temperature. We provide 
an explicit integral expression for (f 2 ) in Appendix |b| 
which depends on A^bs, j(z c ), and the unknown function 
a(z c ). This contribution will simply result in an overall 
scaling of (Ci). A lower bound on this quantity can be 
derived by Jensen's inequality (see Appendix [B]) , which 
yields (f 2 ) > 1/(1 + N ohs a^ ls /6). The behavior of (/ 2 ) 
will generally depend on the sign of a\, since a multi- 
ple collision scenario producing mostly cold spots on the 
CMB will tend to lower To with respect to Tq, driving 
(/ 2 ) > 1 while the opposite is true if a\ > 0. Regard- 
less, we expect (/ 2 ) ~ 1 in the perturbative regime with 
which we are concerned, since N b B ai£i s <C 1 in this case. 

We are now in position to evaluate (C;) in terms of the 
relevant parameters of the theory. Using Eqs. [2l] [32] |4l] 
[44} [45] [46] and performing the sums involving the Gaunt 
coefficients, we arrive at: 



(C ; )=(/ 2 ) 



7r 2 A bsa2^ 2 s / 
21 + 1 {WI 



Ri) + c\ h ( l 



A obs a!^ s A obs (A obs - 



6 



+ 



144 



(47) 



(Not e th at we have dropped the I = piece from 
Eq. (B2) and written out i?°(A D b s ) explicitly as the last 
term?) This is our main result, and characterizes the 
effects of multiple collisions on the CMB power spec- 
trum under our assumptions. The expected angular 
power spectrum depends only on three parameters, A bsi 
a\, and a 2 which encode the model-specific informa- 
tion. The shape of the power spectrum itself results from 
the isotropic and scale-invariant distribution of bubbles 
which should apply regardless of the model potential, as 
long as Hp/Hi 3> 1 and £/ s <C 1. 

As one might expect, the most significant contribution 
to (Ci) is that of wi. From Fig. |6] Ri <C wi so its contri- 
bution is negligible. There is a scale-invariant contribu- 
tion in Eq. |47| proportional to N b S ax, and so one might 
imagine it becoming comparable to that of wi for small 



cti.2- However, this occurs only for very small values of 
«2 such that 



< 



(2/ + 1) 2 
36vr 4 



wi 



(48) 



For example, for the scale invariant piece to be compa- 
rable to that of Wi for I > 10 requires a 2 ^ 10~ n by 
Eq. [48] assuming = .05. For it to dominate at lower I 
would require even smaller a 2 , in which case the effects 
of the collisions would be very difficult to discern except 
for very large N ^ s . Even if it were there, this contribu- 
tion would be degenerate with a slightly different power 
in the primordial perturbations, at least for the 2-point 
function. 

Using this fact, along with (/ 2 ) s» 1, the expected 
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FIG. 7. Expected temperature anisotropy angular power spectrum for multiple collision scenarios from Eq. [49] with iVobs = 100, 
£ is = .05, and for a 2 = 10~~ 6 , 10~ 7 , 10 _s in red, green, and orange, respectively. The Cf 6 are shown by the dotted black line 
for reference. Points in the shaded blue region lie within the approximate WMAP7 error bars [32] ■ The increase in power on 
large angular scales is apparent. 



power spectrum describing a multiple collision scenario 
simplifies to 



(Ci) 



ir 2 N ohs a 2 £,f s 
21 + 1 



wi 



C, 



(49) 



This spectrum is plotted in Fig. [7] for iV bs — 100 and 
various values of a 2 , along with the approximate error 
bars from the WMAP 7 year results [35]. In this scenario, 
the effective power spectrum depends only on a 2 and 
-/V bs- From Eq. 44 the collisions will not contribute any 



statistical anisotropy to the expected spectrum. 

An important feature of the resulting spectrum in 
Fig. [7] is the dramatic increase in power of the low multi- 
poles. This behavior can already be seen in the quantity 
wi as in Fig. [6] and has a physical origin: the flat proba- 
bility distribution over bubble scales yields an equal num- 
ber of large and small bubbles, and the bubble profile is 
"smooth" so that large bubbles contribute only to rela- 
tively low-Z multipoles. Yet larger bubbles both take up a 
larger fraction of the sky (per bubble), and have a larger 
central amplitude. (The larger the affected disk on the 
sky, the farther out on the affected portion of the last 
scattering surface - and hence the closer to the collision 
event - the observer can see; this is evident in the temper- 
ature profile Eq. (21), since the maximum amplitude of 



the temperature perturbation for a given collision scales 
as \T max \ ~ |ai|&,(l - cos^j/2).) Since larger bubbles 
tend to be brighter, searches for collision disks on the 
CMB temperature map, such as those of Ref. [TS], may 
still be able to discern several distinct disks, correspond- 
ing to the brightest collisions, as long as N bs is not too 
large. 

The peak in power occurs regardless of the sign of 
a(z c ), since this effect arises from terms proportional 



to a 2 . For there to be a decrease in power, the scale- 
invariant contribution to Eq. [47] would have to be large 
with evi < 0. Even if this were the case, from Eq. |B14 
the product with (/ 2 ) will still tend to increase (C;) with 
respect to C\ h . 

As iV bs becomes large, the observed spectrum will 
tend to converge to its expectation value Eq. [47] pro- 
viding an increasingly sharp prediction for a given mul- 
tiple collision scenario. This is simply a statement of 
Bernoulli's Theorem. Consider, for example, the one 
point function: neglecting bubbles within bubbles ren- 
ders the individual w\ m statistically independent and a 
multiple collision scenario can then be thought of as nu- 
cleating iV bs bubbles according to the distribution V . In 
this case, Bernoulli's Theorem suggests that with more 
successive nucleation events, the mean value of the vari- 
ous wj m should converge to the expectation value (w^ m ), 
i.e. 



N ohs 



(50) 



as -/V bs — ► co. Multiplying through by iV b s , the sum on 
the LHS of Eq. [50] is simply a particular realization of 
a multiple collision scenario, while the quantity on the 
RHS is the expectation value (wi m ) (c.f. Eq. 37). Thus, 
the expectation value (wi m ) should describe the observed 
wi m with increasing accuracy for large iV bs- Similar ar- 
guments apply to the two-point functions. Therefore, for 
large JV b S ) the observer should expect to see a spectrum 
very close to (C7j). To quantify this assertion and inves- 
tigate the theoretical uncertainty of Eq. [47] across real- 
izations of multiple collision scenarios, one might imag- 
ine computing the variance of the (C;). However, doing 
so would require knowledge of the distribution of a(z c ) 
which we have not attempted to compute and so we do 
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not perform this analysis here. 

Whether or not scenarios with N ODS S> 1 are plausible 
depends on the underlying potential and microphysics. 
From the discussion surrounding Eq. 10 there would 
seem to be cases where a large N oos is quite natural. 
In treating N oos as a free parameter, we have avoided 
the subtleties of this issue and expect the profile Eq. 47 
to provide a good description of the effect of multiple 
collisions on the CMB spectrum for large N ohs . 



C. CMB Constraints On Multiple Collisions: A 
First Look 

Using the results of Eqs. [47j [49j we can preliminarily 
assess how scenarios with N ODS fare in light of observa- 
tions of the CMB temperature spectrum. Of course con- 
straining models based on only one set of observations is 
difficult and a much more in-depth study is required to 
determine how likely it is that we live in a multiple colli- 
sion scenario. We simply wish to answer the question: if 
multiple collisions have impacted our bubble producing 
precisely the expected spectrum, Eq. |49j what do obser- 
vations of the CMB data tell us about the parameters 
iVobs and a 2 ? 

From Fig. [7J we expect that the strongest constraints 
on N b s > 1 in our approximat iorj^l will arise from ob- 
servations of the CMB quadrupole, since its observed 
value is low compared to the best-fit ACDM cosmology 
and multiple collisions will tend to exacerbate this dis- 
crepancy. A detailed statistical analysis in Ref. [33] has 
concluded that the best fit ACDM quadrupole is in fact 
well within the 95% C.L. of the WMAP7 observed value. 
Comparing (C2) with the 95% C.L. upper bound on Cf bs , 
we find that the parameters N a bs and a 2 should satisfy 



N ohs a 2 n c < 3.5 x 10" 



(51) 



Note that this also constrains the magnitude of ol\ , since 
a i < a 2, but not its sign. 

The observed CMB dipole does not constrain the pa- 
rameters N b s , ct2 as strongly as the quadrupole. If we 
imagine that all of the observed dipole has arisen from 
the intrinsic perturbation to the gravitational potential 
from collisions, consistency with the observed value of 
the dipole dictates N b s ct2^c % 1.5 x 10~ 6 . Without 
taking the frame shift into account for \p\ < 1 bubbles, 
the bound would be slightly strengthened, A^ b s Q!2^c ~ 
2.5 x 10 -7 , but still not competitive with that provided 
by the quadrupole. 

A few notes on this result are in order: 



• This result is, of course, dependent upon the de- 
tailed assumptions explained in previous sections 
- namely, the radial profile of a single bubble, the 
flat distribution in bubble angular sizes, the "frame 
shift" , and the approximation that bubbles super- 
pose. However, all of these are unlikely to be dra- 
matically incorrect in a realistic scenario. 

• It is also not necessarily the case that the CMB 
temperature power spectrum is the best observa- 
tional constraint on multiple collision scenarios. 
Observations of e.g. large scale structure, as well 
as CMB polarization and non-Gaussianity, will also 
be affected by the collisions and should be investi- 
gated. These signatures are left for future exten- 
sions of this work. 

• Our results are useful primarily for large N ooa ; for 
a small number of bubbles, searching for disk tem- 
plates on the sky should provide much stronger con- 
straints. 



VI. SUMMARY AND CONCLUSIONS 

In this study, we have begun to address the possibility 
of having many cosmic bubble collisions influencing the 
visible portion of our last scattering surface. This can oc- 
cur despite vanishingly small nucleation rates as long as 
it is compensated for by a large ratio of Hubble constants 
between the false vacuum and the presumed inflationary 
phase in the observation bubble. For such a scenario to 
be compatible with the existence of observers (and to not 
produce effects clearly inconsistent with observation) , the 
effects of the collisions must be mild, suggesting that one 
should be able to treat the effects of collisions perturba- 
tively for models of interest. We were thus led to consider 
a scenario where many bubbles collide with the obser- 
vation bubble while leaving the interior FRW foliation 
intact, supporting inflation and a standard background 
cosmology to the future of the collisions and allowing us 
to study the effects of the collisions on the last scattering 
surface and CMB temperature spectrum. 

Given a number of reasonable simplifying assumptions, 
we are led to a model for these perturbations with the 
following characteristics. 

• The overall effect on the last-scattering surface and 
on the CMB temperature is necessarily small (as- 
sumed for compatibility with observations.) 

• The total perturbation on the last-scattering sur- 
face is a superposition of individual bubble per- 



10 We stress that our analysis applies only to angular scales 1°, 
corresponding roughly to I < 70. To determine the effects of 
collisions on smaller scales requires going beyond the Sachs- Wolfe 
approximation, which we defer to future work. 



turbations (see Sec. IV A) resulting from collisions 
with bubbles decaying to the same vacuum (which 
may or may not match that of the observation bub- 
ble.) 

• Each bubble appears to cover less than half of the 
sky (see Sec. IIIB), with a size distribution that 
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is flat in cos tp, where ip is the angular size of the 
perturbed disk (see Sec. II B I. 



• Each bubble's profile is a 'truncated dipole', with 
central amplitude dependent on ip and z c , a 'kine- 
matic' parameter describing the (dS-invariant) dis- 
tance between the bubbles' nucleation points (see 
Sec.|ilTA|. 



• The details of the kinematics and the inflaton po- 
tential that determine the perturbation strength 
can be folded into two parameters 011,2 ■ The other 
relevant quantity is the total number iV bs of (less 
than full-sky) bubbles on the last-scattering sur- 
face. These three parameters are sufficient to de- 
scribe the effects of many collisions on the CMB 
temperature. 

We put this picture forward as a scenario that is both 
sufficiently simple that predictions can be computed in 
detail from it, and including enough aspects of a realistic 
scenario so that it should give an accurate overall picture 
of the expected effect of a many-bubble scenario. 

The above picture, to a good approximation, yields a 
contribution to the / < 70 power spectrum of 



(Ci, 



ir 2 N ohs a 2 S,l 
21 + 1 



-wi, 



(52) 



where £; s ~ 2^0^, a 2 encodes features of the inflaton po- 
tential that determine the amplitude of a given collision 
perturbation, and wi is a fixed function given by Eq. 42 
and shown in Fig. [6] This equation holds with increas- 
ing accuracy as iV bs becomes large. The resulting power 
spectrum contribution is very red-tilted (see Fig. [7|) , so 
that the strongest limit on a multi-bubble scenario comes 
from the quadrupole. 

In our universe, a multiple collision scenario with iV bs 
visible bubbles and average brightness described by a 2 
consistent with observations of the CMB power spectrum 
would require roughly a 2 il c ^ 3.5 x 10~ 9 /N ^ s . This 
bound can in turn be used to constrain a given well- 
specified scenario by developing in detail the link between 
the inflaton potential and the parameters a 2 , and 
N ohs . 

We emphasize that this is a first study aimed at laying 
the groundwork for such scenarios and computing con- 
straints using the power spectrum. It is quite possible, 
and even likely, that other observables like polarization 
or non-Gaussianity could give significantly stronger con- 
straints, or provide a cleaner method for the actual de- 
tection of bubble effects. This set of methods comprises 
a complement for those already developed to investigate 
scenarios with a few bubble impacts. In either case, re- 
finement of these tests will soon provide real constraints 
on the parameter space of eternal inflation, and might 
even produce evidence for eternal inflation itself. 
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Appendix A: Intervening Bubbles 

In this Appendix we detail the calculation of the ex- 
pected number of intervening bubbles to intersect a col- 
liding bubble with nucleation center (z c , p). Since we are 
interested in intervening bubbles intercepting late-time 
collision bubbles, we can neglect the angular dependence 
of the initial value surface and choose 9 C = 4> c — with- 
out loss of generality. 

By the discussion in Sec. |IV A[ we must calculate the 
4-volume I(z c ,p) with the cosmologies fixed. T(z c , p) is 
given by 



T(z c , p) =Vf LC [r] (z c ,p),T (z c ,p)} 

- Vf LC [ Vc (z c ,p),T c (z c ,p)} - Vf LC [0,0] 

(Al) 

where Vf LC [r], T] denotes the 4-volume in the past 
light cone (PLC) from the point (77, T) and (770, To) is 
the point of intersection of the future light cone from 
(rj c (z c , p), T c (z c , p)) and the observation bubble wall. 

It remains to determine whether a point x 2 — 
{ r l2,T 2l 6 2 ,4) 2 ) is within the past light cone of x\ — 
(t}\,T\,9i,4>i). In the embedding space, two points X^, 



X 2 are timelike or null separated if and only if 



(A2) 



Restricting this relation to the dS hyperboloid gives an 
equivalent relation 

sin 771 sin r/ 2 

P(xi,x 2 ) = - tanTi tanT 2 + _ ^ m u) i (xi)Uj(X2)d i j 



cos Ti cos T 2 



cos 771 cos r) 2 
cos Ti cos T 2 



> 1 



(A3) 

for timelike or null separated points, where the u>i are 
defined as in Eq. [I] We can then numerically integrate 
over the 4-volume to the past of each of the points Xi in 
Eq. Al weighted by the step function Q(P(xi, x) — 1) to 
obtain T(z c , p). The integration within the set of points 
PLC[0,0] can be done analytically, however the angular- 
dependence of the lightcones away from the origin makes 
it difficult to do so for the other points. We cut off the nu- 
merical integration at T corresponding to the intersection 
of the PLC from (ry c ,T c ) and that from the origin. This 
does not count regions near past infinity (corresponding 
to early time bubbles), however these portions of the false 
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vacuum are typically included in PLC[r] c , T c ] <1PLC[0, 0] 
and so will not contribute to A^„ t . 

To obtain A; nt we use Eq. [20] with a particular choice 
of A. In Fig. [5] we use A such that iV bs = 1 from Eq. [10] 
Increasing iV b s will correspond to increasing either A or 
£; s , however we still find Ni n t <C 1 for cases of interest. 



metric about the z-axis, only the ml = modes con- 
tribute and we are left with 



where 



, w i'0' D 0'm (<Xi,Pi,nti)Sw 



(B2) 



Appendix B: Calculating Expectation Values Of 
Correlation Functions 



In this Appendix we detail the calculations of the ex- 
pected one- and two-point correlation functions Eqs. 41 



44 Our conventions follow that of Ref. [30 a , to which 
we refer the reader for further details as well as useful 
identities used in the calculations we present. 

We are interested in the two-point function of the tem- 
perature anisotropy which in turn will involve the one- 
and two-point functions of uii m . Our analysis is simpli- 
fied by the following observation: for a given function g 
expanded in spherical harmonics in a coordinate system 
(primed indices), the coefficients in a rotated coordinate 
system (unprimed indices) are given by [2^1 130] 



9h 



91' 



(Bl) 



where 7, 9, <fi are the Euler angles specifying the transfor- 
mation between the two coordinate systems and D m , m is 
the Wigner rotation matrix. Since wi m is a sum of the 
coefficients corresponding to the individual bubbles, we 



can write wi r 



+N 2 *>l. 



and evaluate each w\ m in 



a frame in which the collision is centered at the north 
pole 9 = 4> = and rotate them back into the original 
coordinate system via the transformation Eq. |B1| 

Since, in the polar (primed) frame, the profile for each 
bubble on the sky by construction is azimuthally sym- 



I'm' 



y; m ,(MK(fi)k=o,0 Ol =o^ (B3) 



D 0'm ("t, A)7t) 



■iy 



Air 
21 + 1 



YL (A, 7.) (B4) 



In our calculations it will often be convenient to con- 
sider the quantity wi>o>(ip), corresponding to the coef- 
ficient wi m from a collision centered at the north pole 
9 = 4> = (hence the prime), with angular scale ifi, and 
with the strength a£i s factored out: 



1 

27ra<£j. 



J l'0'- 



(B5) 



This is precisely the quantity explicitly defined in Eq. 43 
The relevant Euler angles to rotate from the polar frame 
to the original frame are given in our convention by = 

9oi, 7» = 7T - (f> i (A = 7T - 0Q», 7« = n + 0Oi) for * < ^1 

(i > N\). The coefficients ai m are related to wi m via 
dim = fwim where we have neglected the monopole piece 
arising from the I = integral. 

To carry out the ensemble averages for the various cor- 
relation functions, we average by the density V defined 
in Eq. |36] We use the shorthand 



i . v 

sm — 

16tt 2 



(B6) 



for brevity. In this prescription, the expected one point 
function is given by 



H)" 1 ^ <<*(*,)) [Ni (Y lm (9, tt - 0)) (w m + N 2 (Y^n - 

16?r 3 



,ir + (f>)) (wi> >(-ip))^J S w 



^V"obB«l^a(-l) 



21 + 1 



Sw / <w(V0ff(V0# / Y lm (9,-K-^)dn 



v / 2T+T 



Sii'SioSmo / wii >(ip)g(ip)dip 



(B7) 



Here the brackets ()^,< 7r (> 7r ) denote averaging over %p in 
the interval [0, tt] ([tt, 2tt]). In going from the first to sec- 
ond line, we use Eqs. |B2[|B4| and |B5| to express the sum 
as a collection of averages in the polar frame with the ap- 
propriate rotation coefficients. In going from the second 



to the third line we use the statistical independence of the 
collisions write out the averages explicitly in terms of the 
averages for an individual collision (as per Eq. 37 1 and use 
the fact that / Y t * m (9, tt - ^dil = J Y{ m {jt - 9, tt + (j))dfl 
to combine the Ni and N 2 terms. The last line follows 
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from integrating the spherical harmonic. 

Only the monopole part contributes to the one-point 
function and we can perform the integrals over ip, using 
Eq. (43) to express wq>q/ as 



Inserting this into the last line of B7 and integrating 
yields Eq. (Ell). 



™0'0'WO 



1 



sin 4 -6 (tt - ib) + cos 4 -6 (ib - tt) 
4 ' 4 



(B8) 



Moving on to the expectation value of the two point 
function, (^^ui^J, we compute: 



{wi imi w hm2 ) = (-1) 



y/(2k + 1)(2Z 2 + 1) v , 

vmrmth fc {a{Zc)2) h^k^) < y ^*w 



+ ^2(a(z c )) 2 (w^viipdw^viipj)) (*£ mi (0 i( tt - <bi)Y l2rn2 {6 h ir ~ fa)) 

16n 3 £? s 5 h i> 5 h i> i f 2 

y/ {Ml + lj(z« 2 + 1J v J 
+ N ohs {N obs - l)a\ (w l[0 ,) (wi> a0 >) (yjtmx) 0-W 

21 



(B9) 



+ 1 ^l^libVmi +^obs(^obs - 1) I g I 



I 

In the first line, we write out the average as a collection ensemble average. Let us first consider the quantity / _1 
of iV bs collisions in the polar frame with the appropri- and write 



sky • 



(BIO) 



ate rotation coefficients. Note that the function wi>o'(ij>) 
is real-valued, hence only the complex conjugate of the 
spherical harmonic appears here. The second line follows 
from the commutation of the averages. In the third equal- 
ity, we use the statistical independence of the individual Th e sky average for a given Wi is independent of the 
collisions to write out the average two-point function ex- collision's location on the sky and thus depends only on 
plicitly in terms of the averages for an individual collision We can thus define the a-independent quantity 

and integrate over a(z c ). We don't write out the one- ^ . 1 

point averages explicitly for brevity; they are precisely S(ip) = / u/(n) 

as in Eq. |B7| Between the third and fourth equality we 
perform the remaining integrals, again combining the N± 
and N 2 terms, since (YJ* mi (0,7T- faYi 2m2 (9,-K - fa) = 
( Y i* imi (K ~ 0,n + (b)Yi 2 , n2 (n - 8,n + fa) and similarly 
for the one-point averages. 



in terms of which, f 2 is given by 
f= (l + 



2,fn 



wo'o'O) 
(BH) 



(B12) 



Finally, we can compute (f 2 )- The quantity / relates 
the average CMB temperature to the average tempera- Using these defi niti ons, squaring, and averaging by the 



ture without any collisions. One must then average the distribution Eq. B6 we arrive at our expression for (f 2 ) 
spatial temperature profile over the sky, then take the in terms of iV bsj ot[z c ), and £; s : 



/,2\ _ f (^Tr) 2 g(^i)9{ip2)j(z cl )j(z C 2)dTpidip 2 dz cl dz C 2 

U ; ~ J 1 + N Uu (a(^)5(^i) + a(z c2 )S(ip 2 )) + AUs(a(*ci)6s) 2 <S 2 (>i) + iW^obs - l)a(z cl )a(2 c2 )£2 <S(^i)<S(V> 2 ) ' 

(B13) 



A useful lower bound on (f 2 ) can be derived from Jensen's inequality, whereby (f 2 ) > (1//) 2 ■ Dropping 
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the term quadratic in Nai£i s results in 

(/ 2 >>(l+ iVob f 1 ^)" 1 . (B14) 



In our perturbative approach we expect (/ 2 ) ~ 1. This 
quantity simply results in an overall scaling of the power 
spectrum. 
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